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A developing surface frontal system with an overlying polar front 
jet and subtropical jet stream merging region was present over the con-=- 
tinental United States during 22 to 24 March, 1966, Possible cross=- 
over flow between these two jet streams has been studied in detail by 
isentropic air trajectories at four potential temperature levels from 
31@ K to 335 i: 

The results of the trajectory analyses on the four 6 surfaces 
have indicated a cross-over of air parcels from the subtropical jet 
stream to the polar front jet stream, No evidence has been found for a 
reciprocal relationship for air parcels initially associated with the 


polar front jet between 310 K and 335 K isentropic levels, 
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SECTION I 


INTRODUCTION 


Studies by Reiter and Whitney (1965) of the region where a polar 
front jet and a subtropical jet stream are in close proximity to each 
other have produced a provecative hypothesis that there exists a cross- 
over flow between the two respective jets, One corollary is that cer- 
tain air parcels in the polar front jet, rather than follow a normal 
cyclonic trajectory over the continental United States, crosseover and 
proceed along a trajectory path similar to that of air parcels in the 
subtropical jet stream, A second corollary proposes that, rather than 
follow a basically anticyclonic trajectory, certain air parcels in the 
subtropical jet crosseover and then continue along a trajectory path 
similar to that of air parcels in the polar front jet stream flow, 

Regardless of the validity of the hypothesis, the isentropic tra- 
jectory of a representative air parcel in either jet can be computed by 
the method developed by Danielsen (1961), This computational method 
describes with considerable accuracy the movement of air parcels along 
an isentropic surface, However, unlike several other approaches to 
calculating air trajectories which assume that air parcels move quasi- 
horizontally on a surface of constant height or pressure, the Danielsen 
approach gives the necessary third dimension -— a measurement of the 
average vertical motion, 

Throughout this study diabatic effects on trajectory motion in 
this non=precipating region are neglected because diabatic contribu- 
tions to vertical motion are normally small compared to adiabatic con-= 


tributions, 
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During late March 1966 a strong developing surface low with as- 
sociated frontal systems moved east=northeast across the continental 
United States, Also associated with this surface system were a strong 
polar front and a subtropical jet stream in close proximity to each 
other as fi teat rated i Figure 1, The period between 22 March 0000 GMT 
to 24 March 0000 GMT was selected as being most favorable for testing 


the hypothesis of Reiter and Whitney (1965), 
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SECTION I1 


NATURE AND ACQUISITION OF THE DATA 


All the data in the analysis comes from radiosonde and rawinsonde 
checked data found in the U.S, Weather Bureau Northern Hemisphere Data 


Tabulations for 22, 23, and 24 March 1966, 





The period under study contains five observation times, 22 0000 
GMT, 22 1200 GMT, 23 0000 GMT, 23 1200 GMT, cat 24 0000 GMT, For each 
observation time, an area of interest is determined by plotting an area 
encompassing not only the surface system but also the region ot the 
merging jets, Table 1 shows the data areas for each observation time, 
Having delineated the regions of interest, the checked data is utilized 
from each upper-air station within these regions, 

The data for each upper-air station consists of all standard pres- 
sure levels in 50 mb intervals from 1000 mb to 200 mb, From 200 mb to 
100 mb the interval is 25 mb, Temperature in degrees Celsius, heights 
in meters, Wind direction in degrees, and wind speed in meters per sec- 
ond are tabulated for each standard level, In addition to the standard 
levels, the significant levels are recorded with their tabulated temp- 


eratures and heights 
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Table 1, 


Data areas used for each observation time, 


OBSERVATION | 


TIME 





22 March 
0000 GMT 


22 March 
1200 GMT 


23 March 
OOOO GMT 


23 March 
1200 GMT 


24, March 
0000 GMT 


U,S,-Mexico 
Border and 
Gulf Coast 


U,S,—Mexico 
Border and 
Gulf Coast 


U,S,—Mexico 
Border and 
Gulf Coast 


U,S,—Mexico 
Border and 
Gulf Coast 


U,S,—Mexico 
Border and 
Gulf Coast 


MW 


West Coast 
U.S, 
115°W 
105°w 


105°W 


105°W 


BAST 


90°W 


80°W 


75°W 


75°W 


65°W 


SECTION IIl 


THE MONTGOMERY STREAM FUNCTION COMPUTER PROGRAM 


=O 
Assuming de/dt , the equations of motion in © coordinates become 


gu 4 udu du -_9M ify 
9 ‘ aXe pe 4o aX» (1) 
av +ugv yyau - 9 =afh , (2) 


where the isentropic stream function is defined as 
M=cpl+9z., (3) 


For the air trajectory method as derived by Danielsen (1961) and 
later outlined by Mahlman (1965), the Montgomery stream function is 


given by 


Vive =Cple + 92p 1 Polen Qn B) ° (4) 


As observed by Danielsen (1959), large errors of M result if val- 
ues of I, and Po are selected independently of each other from the 
plotted soundings, Consequently, for this program, the dependent re-= 
lationship of T, and Pg on an isentropic surface, as given in Poisson’s 


equation, allows M to be written 


M = Cp le + ql ts RT hf Be (2) yn) ; (5) 


as shown in Mahlman and Kamm (1965) M is now in a convenient compu- 
tational form which can be handled readily by the computer, 


The basic program is presented in Appendix I and will be briefly 
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outlined in this section, After reading in the data, all temperatures 
are converted to degrees Kelvin, Potential temperatures are next com- 
puted for each level, For this study, a range of potential temperatures 
between 290 K and 350 K at 5° intervals is assigned for computation, 

The data or ae searched for the first @ level above and below the 
predetermined 6 level, With these two 6 values and the predetermined 
@ value, T, is computed, Pp is then computed from Poisson’s equation 
using Tp and ©, Finding next the first standard levels above and below 
the arbitrary isentropic surface and calculating T by summing the weight- 
ed mean temperatures between all data points from the lower standard 
level to the arbitrary @ level, the Montgomery stream function is com 
puted in three segments according to equation (5), Finally, the wind 
direction and wind speed, # interpolated according to the procedure used 
by Mahlman and Kamm (1965), After completing the calculations for one 
particular isentropic surface at all stations in the data area, the pro=— 
cedure is repeated for the next higher relevant @ level, 

The final computer output gives for each isentropic surface the M 
value, pressure, static stability, wind direction, and wind speed for 
all reporting stations, With these five parameters, the isentropic 


trajectories can be computed and checked for accuracy, thus enabling 


the three-dimensional tracking of selected air parcels, 


iy 


SECTION IV 


DETERMINATION OF RELEVANT ISENTROPIC SURFACES 


Possessing now the computed parameters M, Po, wind direction and 
wind speed for all the isentropic surfaces at 5° intervals between 290 K 
and 350 K, the next step in the procedure is to determine particular 6 
levels which are representative of the isentropic surfaces present in 
the polar front jet and the subtropical jet streams, 

The observation time 23 March 12000 GMT is selected as the synoptic 
map most representative of the five observation times, A north-south 
cross section on longitude 95° Wis constructed with plotted values 
of wind speed and potential temperature as depicted in Figure 2, An- 
alysis of this cross section shows that four isentropic surfaces might 
possibly be of interest —— 310 K, 320 K, 330 K, and 335 K, These four 
6 surfaces are chosen because of their physical location with respect 
to both jet streams and their heights above sea level, so that there 
will be a minimum of orographic and lower level interference, 

Basic maps are plotted on 1:15,000,000 scale charts for the four 
@ surfaces, Hach of these isentropic surfaces is analysed for the 
Montgomery stream function using a 0,5x10/ cm” sec™* interval, for pres-= 
sure using a 50 mb interval, and for wind speed using a 10 m sec! iso- 


tach interval, Figures 3, 4, and 5 are typical examples of these anal- 


yses, 
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SECTION V 


COMPUTATION OF ISENTROPIC AIR TRAJECTORIES 


Mahlman’s outline (1965) of the Danielsen method (1961) develops 
in detail the two working equations for computation of isentropic air 


trajectories, 

t, rg wy 
Japa Oe Md rE (VeVi), (6) 
7 


D=\VAt. (7) 


ol 


The left hand term in equation (6) represents the local change in 
the stream function, The local change is graphically determined by an- 
alysing the 12 hour local changes of oM at each upper-air station or 
data point, For the four time periods of this study, aM is analysed 
using an increment of 0,25x10? cm* sec™*, oM is then averaged along the 
trajectory path to obtain an approximation to the left hand side of 
equation (6), Figure 6 illustrates the local change of M for the period 
22 March 0000 GMT to 22 March 1200 GMT on the 330 K 6 surface, 

The second term on the right hand side of equation (6) expresses 
the change in kinetic energy of the horizontal wind, After converting 
the units of this kinetic energy term to the units of M (107 om” sec”) 


and resarranging terms, equation (6) can be written as 
2 é 
M,=M, + 0.0008 CV-V,),+AM (8) 


pi 


for V expressed in m sec” 
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The units of equation (7), the distance equation, are converted to 


degrees latitude for a 6 hour trajectory, becoming 
D=0.19 Vy « (9) 


The objective now is to locate a value of M at the final point on 
a trajectory path such that equations (8) and (9) are simultaneously 
satisfied, Detailed procedures for calculating isentropic air traject- 
ories can be found in Danielsen’s paper (1961), However, the procedures 
used for this study will be briefly outlined for one time period, 

For the first 12 hour period, 22 0000 GMT to 22 1200 GMT, a rep- 
resentative sample of initial data points in and around the isotach 
center of the polar front jet and the subtropical jet streams is selected 
on the 22 0000 GMT M analysis for the four relevant 9 surfaces, The 
first estimate of the path of the trajectory is obtained by sketching 
both the 22 0000 GMT and the 22 1200 GMT streamlines around one of the 
initial data points, An average streamline is determined by giving 
more weight to the 22 0000 GMT streamlines for the first 6 hour period 
and to the 22 1200 GMT for the second 6 hour period, A 6 hour traject-— 
ory is computed by equation (9) using the wind speed at the initial data 
point and is plotted along the average streamline using i increments, 

A mean wind at the 6 hour position is determined from the 22 0000 GMT 
and 22 1200 GMT isotach analyses, This mean wind is then averaged 
with the initial wind to determine an estimated mean wind over the 
first 6 hour period, Equation (9) is again employed to calculate a 
better eStimate of the 6 hour trajectory position, 


Using the mean wind at the 6 hour position, the trajectory for the 
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next 6 hour period is computed by equation (9), The winds at the 12 
hour and 6 hour positions are averaged and this resulting mean wind is 
used to re-estimate the 12 hour position, 

Now, the trajectory computed by using equation (9) must be tested 
by equation (8), To, calculate OM, the trajectory determined by the 
distance equation is superimposed on the respective aM analysis and 
divided into four equal segments, Simpson’s rule is utilized to calcu- 
late the 12 hour local change of M integrated along the trajectory path, 
Having aM, V1; Vo, and M,, M, is computed, 

At this stage in the procedure, an observed value of Mo, determined 
from the first estimate to the trajectory and the 22 1200 GMT M analysis, 
and a computed value of Mj, determined by equation (8), have been ob- 
tained, M>_opbs is then compared to Mo comps If their absolute difference 
is equal to or less than O,1«10? cm? sec”, the initial trajectory sat- 
isfies both equations and this trajectory represents the 12 hour air 
parcel movement on the isentropic surface, 

Generally, for the first estimated trajectory, | Mo comp = 5, 
is greater than 0,1x10’ cm? sec”, If this is the case, the 12 hour 
position calculated by equation (9) is moved perpendicular to the mean 
flow towards the M5 comp value obtained from equation (8), A new V5 at 
this second estimated 12 hour position is determined from the 22 1200 
GMT isotach analysis, 

The initial trajectory is now adjusted throughout its entire path 
to conform to the second estimated 12 hour position, The resulting 
differences between the adjusted trajectory and the initial trajectory 
have changed the 6 hour position, Hence, a new mean wind over the first 
6 hour period and the second 6 hour period is calculated by the method 


outlined previously, 
AA 


Following the original procedures with the new 6 hour mean winds, 
a second estimated 12 hour trajectory position on the adjusted trajec- 
tory is computed, If — - Mo obs | new 28 less than or equal to 
O,1x10/ om* sec", the second estimated trajectory satisfies the dis- 
tance and the energy equations simultaneously, If not, the previous 
procedures are again repeated until convergence is reached, Thus, an 
air parcel has been successfully tracked for 12 hours on an isentropic 
surface by successive approximations, 

Figure 7 illustrates the final estimated trajectories for the 


22 0000 GMT to 23 OQ000 GT time period on the 330 K isentropic surface, 
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SECTION V1 


RESULTS OF THE ISZNTROPIC AIR TRAJECTORY COMPUTATIONS 


After computing a representative sample of trajectories for each 
of the four isentropic surfaces, (310 K, 320 K, 330 K, and 335 K), for 
all four time periods, the trajectories are superimposed on the respec— 
tive pressure analyses, so that Po can be determined for the initial and 


final position, ® is then calculated for each 12 hour trajectory, where 


te 
wat fae dt = CFL Py nbz nes . (20) 
' 

To illustrate the three-dimensional air parcel movement, an analy- 
sis of trajectory position relative to the surface front position is 
made, The positional relationship between a trajectory and the moving 
surface frontal system, as shown in Figure 7, appears quite consistent 
for all four 6 surfaces, Trajectory 2 on the 330 K surface, initially 
over Las Vegas, Nevada (Station 386), has an average downward vertical 
motion (4 positive) of +16 mb/l2hrs for the 22 0000 GMT to 22 1200 GT 
period, This downward vertical motion is reasonable since the Trajec- 
tory 2 position is generally lagging the cold front position for the 
12 hour period -——- thus remaining in a subsidence region, 

As Trajectory 2 progresses to the northeast, its position relative 
to the frontal system is, on the average, leading the surface system, 
This leading situation is indicative of the average upward vertical 
motion of -1 mb/l2hrs for the 22 1200 GMT to 23 0000 GMT period, 

This leading and lagging relative relationship between isentropic 
air trajectories and the surface frontal system was consistent for all 


time periods and isentropic surfaces, Thus, through the laborious 4 
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level computations, a physical understanding of air trajectory three- 
dimensional motion has been attained, In addition, a mental picture of 
this motion can be generated by relating these trajectories to the sur= 


face phenomenon, 
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SECTION VII 


TESTING THE HYPOTHESIS 


Recent studies by Danielsen (1961), Reiter and Whitney (1965), and 
Mahlman (1965) have investigated cases of apparent jet stream cross-over 
at isentropic levels below 310 K, Assuming the crosssover phenomenon 
does exist at the lower levels of the atmosphere, it is reasonable to 
inquire as to the existence of a compensatory crosseover flow above 
the 310 K isentropic level, Therefore, this hypothesis will be tested 
by further analyses of trajectories associated with the respective jet 
streams within the merging region, 

The isentropic air trajectories are analysed for their positions 
relative to isotach centers of the polar front jet and the subtropical 
jet streams, The trajectories are further analysed for their positions 
relative to the region of apparent jet stream merging, As illustrated 
in Figure 38, Trajectory 1 at 22 0000 GMT is located within the 50 m 
sect isotach center of the polar front jet stream on the 335 K sur- 
face, Maintaining a cyclonic curvature with an average downward vertical 
motion of +8 mb/l2hrs for the succeeding 12 hour period, an air parcel 
associated with Trajectory 1 remains in the polar jet flow and by 22 
1200 GMT is well out of the merging region, Referring again to Figure 8, 
an air parcel associated with Trajectory 2 follows a similar path as 
that of Trajectory 1, but is displaced more to the south throughout 
the 22 0000 GMT to 22 1200 GMT period, As before, this trajectory can 
be tracked through the merging region and on to the northeast, 

After analysing the trajectories associated with the polar front 
jet stream on the four 9 surfaces for the four time periods, the mean 


three=dimensional flow of air parcels in the polar front jet is deter- 
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mined to be basically cyclonic in nature -—=- descending when moving south- 
east and aScending when moving northeast, This mean flow corroborates 
previous studies, 

For the 22 0000 GMT to 22 1200 GMT period, the apparent jet merging 
region is over New Mexico and Arizona, Trajectory 3 of Figure 8 at 
22 0000 GMT begins within the merging region over Tucson, Arizona 
(Station 274), An air parcel following this trajectory appears to 
disassociate from the main flow of the subtropical jet and proceeds to 
the northeast with an average upward vertical motion of -10 mb/]2hrs 
for the first 12 hour period, By 22 1200 GMT, the air parcel associated 
with Trajectory 3 has entered the isotach center over Missouri and 
continues to move northeastward toward the Great Lakes region, This 
phenomenon of apparent disassociation of air parcels from the subtropical 
stream flow is evident on the four 6 surfaces for the trajectories with 
initial position over Tuscon at 22 0000 GMT, 

Trajectory of Figure 8 is located over El Paso, Texas (Station 
270) at 22 0000 GMT, Unlike Trajectory 3, an air parcel associated with 
Trajectory 4 remains within the anticyclonically curving subtropical 
jet stream flow for the succeeding 2/4 hour time period, 

To substantiate the statement made in the previous paragraph con- 
cerning apparent disassociation of air parcels from the subtropical jet 
stream flow on a © surface other than 335 K, refer again to Figure 6, 
Trajectory 3 on the 330 K surface also begins over Tucson, Arizona at 
22 0000 GMT, This trajectory path appears to crosssover and by 22 1200 
GIS possesses the characteristic flow pattern of the polar front jet 
stream, 

By 23 March 0000 GMT, the jet stream merging region is located over 


Texas, Trajectory 2 on the 330 K isentropic surface in Figure 9, 
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initially located over Del Rio, Texas (Station 261), appears to continue 
within the subtropical jet influence and maintains an average upward 
vertical motion of -29 mb/l2hrs for the 23 0000 GMT to 23 1200 GMT 
time period, 

The isotach analyses of the four isentropic surfaces for the 23 
1200 GMT to 24 0000 GMT period defines the merging region of the polar 
front and the subtropical jet streams most clearly, Trajectory 3 of 
Figure 10 was selected because not only was the air parcel following 
this trajectory path associated with the subtropical jet flow, but the 
23 1200 GMT position over Midland, Texas (Station 265) was well within 
the merging region, Trajectory 3 on the 320 K and 330 K isentropic 
surfaces appears to crossover and proceed cyclonically to the north- 
east, apparently identifying itself with the polar front jet flow at 
these levels as illustrated in Figure 10, However, trajectories on the 
335 K and 310 K surfaces located over Midland at 23 1200 GMT are dis- 
placed to the south of the 320 K and 330 K trajectory position, indicating 
a continued association with the subtropical jet stream flow, It is 
interesting to note that Trajectory 4, located initially on the anti- 
cyclonic side of the subtropical jet stream, apparently remains well 
within this flow, 

The subtropical jet is located north of its normal position for the 
23 March 1200 GMT to 24 March 0000 GMT period, Trajectory 3 of Figure 
10 is initially located on the cyclonic side of the jet and in the 
central section of the merging region, Trajectories computed to the 
soutn, or on the anticyclonic side of the subtropical jet stream, follow 
the classical anticyclonic path over the southern United States as 


shown in Figures 7 and 9, 


BD 


Potential vorticity (P) of an air parcel is known to be a conserva- 
tive property for adiabatic and frictionless air motions, As a further 
test for the four trajectories with apparent crosseover flow, P is com 
puted for the initial, 12 hour, and 24 hour isentropic trajectory posi- 
tions, If the three P values do not vary appreciably, this lends further 
credence to the validity of the trajectories indicating cross#over flow, 
Consequently, the fact that potential vorticity is conserved, as shown 
in Table 2, helps substantiate the results of the isentropic trajectory 
computations, | 

As shown in Figure 8, Trajectory 3 on the 335 K isentropic surface ~ 
for the first 12 hours follows a cyclonic path and for the second 12 
hours an anticyclonic path, If the assumption of conservation of po- 
tential vorticity is valid for this assumed frictionless and adiabatic 
motion, then the large 24 hour relative vorticity change which reflects 
this curvature change must be explained by a significant increase of 
one of the parameters in the P equation, One possibility is a rapid 
increase in the Coriolis parameter, but for this trajectory, variations 
in f are small as shown in Table 2, However, Table 2 shows that the 
relative vorticity decrease is compensated by a rapid static stability 
increase, To some extent this argument helps justify the physical 


mechanisms involved in the crosseover phenomenon, 
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Table 2, Approximate values of potential vorticity computed at rela- 
vant trajectory positions, Note the static stability and absolute vor- 
ticity values corresponding to respective trajectory curvatures shown 
in Figures 7, &, and 10, 
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SECTION VIII 


CONCLUSIONS 


As indicated in Figures 7, 8, and 10 and their subsequent analyses, 
crosssover of air parcels in the subtropical jet stream flow to the 
polar front jet flow apparently exists between the 310 K and 335 K 
isentropic surfaces, This conclusion is dependent upon the validity of 
the adiabatic assumption, However,refercing again to Figures 7, 8 
and 10, it can be observed that small directional wind shears in the 
vertical exist through the 10° isentropic layer, Therefore, diabatic 
effects, even if present, would not appreciably affect the results, An 
additional implication from these analyses is the existence of a mass 
transport phenomencn in the jet stream region as observed by Reiter and 
Whitney (1965). 

The cross#over phenomenon of air parcels initially in the polar 
front jet stream is not observed in this case study for the relevant 
isentropic surfaces, However, the apparent lack of air parcel cross- 
over to the subtropical jet stream flow at these © levels might be ex- 
plained by the limited number of isentropic surfaces, synoptic situa- 


tions, and parameters utilized in this research, 


Al 
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APPENDIX I 


The Montgomery stream function computer program developed for this 


study (Sections A and B), 
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APPENDIX I; Section A, List of symbols used in the computer program, 


SYMBOL MEANING 

DD1 wind direction at first standard level below TH1 
DD2 wind direction at first standard level above TH1 
DDTH interpolated wind direction at TH1 

F Montgomery stream function 

Pte first segment of F 

2 second segment of F 

F3 third segment of F 

FFL Wind speed at first standard level below TH1 
FF2 Wind speed at first standard level above TH1 
FF TH interpolated wind speed at TH1 

PP pressure at first standard level below TH1 
PrZ pressure at first standard level above TH1 
PPS pressure at first level below TH1 

PEE pressure at first level above TH1 

PTH pressure at TH1 

STAB stability parameter = - sS 

T temperature 

7B T at first level below TH1l 

TBAR average T of layer between PP and PP2 

TH potential temperature 

THB TH at first level below TH1 

THe assigned TH value 

Tea TH at first level above TH1 

Ti T at first level above TH1 


AA 


SYMBOL MEANING 
TTH T at THL 


ZTH geopotential height of TH1 
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The Montgomery stream function computer program, 


Section B, 


APPENDIX I; 
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SEARCH DATA FOR VALUES ABOVE AND BELOW LEVEL BEING COMPUTED 
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COMPUTE PRESSURE AND TEMPERATURE AT LEVEL OF INTEREST 


TTH TB + (TT - TB) * C( THI - THB) / (THT - THB)) 
PTH 1000.0 * ({(TTH / TH1L) ** 3.5001) 


FIND NEAREST STANDARD LEVEL 


IPP = PTH / 50.0 

AAA 2 IPP * 50 

fF (G(s) ~EQ. AAA) GO TO 61 
— + 

GO TO 60 

pp = Gttyl) + 50.0 

TF (GJo1) .EO. PP) GO TO 63 

G0 T@ 62 

PP = GlJel) 

HGT = GlJy3) 

FFl = G(J,5) 

DD1 = G(J,4) 

AA = 0.0 

IF (PTH —- G(J+#1,1)) Tl, T2>9 T2 

AA = AA ‘ (TUS) + TOS41)) * (GO0,1L) - Gld4],1)) 
= + 

G0 TO 70 
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COMPUTE MONTGOMERY STREAM FUNCTION IN THREE SEGMENTS 
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SO OOM AGO 


COMPUTE MONTGOMERY STREAM FUNCTION 
Fo = Fl + F2 + F3 

COMPUTE HEIGHT OF THETA SURFACE 
ZTH = HGT - (2.8704 * TBAR / 9.8) * (ALOG(PTH) 
650) Iy STABy Fl, F2y F3y F 
Ky {10s 5X2 6032 SX, 3P4EL 


X 5. 
THI PTH, TtH, 2TH, DOTH 
J 


5) 
Xo 6F10eL, ///) , 


NNOZOONLNS 
zA0O O00R702X 
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